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Quantum discord, as a measure of all quantum correlations, has been proposed as the key 
resource in certain quantum communication tasks and quantum computational models without 
containing much entanglement. Dakic, Vedral, and Brukner [Phys. Rev. Lett. 105, 190502 (2010)] 
introduced a geometric measure of quantum discord (GMQD) and derived an explicit formula for 
any two-qubit state. Luo and Fu [Phys. Rev. A 82, 034302 (2010)] introduced another form of 
GMQD and derived an explicit formula for arbitrary state in a bipartite quantum system. However, 
the explicit analytical expression for any bipartite system was not given. In this work, we give out 
the explicit analytical expressions of the GMQD for a two-parameter class of states in a qubit-qutrit 
system and study its dynamics for the states under various dissipative channels in the first time. 
Our results show that all these dynamic evolutions do not lead to a sudden vanishing of GMQD. 
^ " Quantum correlations vanish at an asymptotic time for local or multi-local dephasing, phase-flip, 

| and depolarizing noise channels. However, it does not disappear even though t — > oo for local 

trit-flip and local trit-phase-flip channels. Our results maybe provide some important information 
for the application of GMQD in hybrid qubit-qutrit systems in quantum information. 
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^"1 1. Introduction 

Quantum systems with quantum correlations, have some fundamental applications for quantum information process- 
1 ing. However, entanglement does not describe all the aspects of the quantum correlations exhibited by a multipartite 
quantum system. There are some other kinds of quantum correlations without entanglement that are also responsible 
for the quantum advantages over their classical counterparts [H4To| . Aiming at capturing such correlations contained 
in a bipartite quantum state, Ollivier and Zurek (Til ] introduced quantum discord as a measure of all quantum corre- 
. ■ lations in 2001. Recently, quantum discord attracted a lot of attention. Unfortunately, quantum discord is difficult 
to calculate and it is usually based on a numerical maximization procedure. There are few analytical expressions, 
including some special cases G^HHj]. For general two-qubit mixed states, the situation is more complicated. To avoid 
this difficulty, Dakic et al. [l7| ] introduced a geometric measure of quantum discord (GMQD) which measures the 
quantum correlations of a quantum system in a given state through the minimal Hilbert-Schmidt distance between the 
given state and a zero-discord state, and they derived an explicit formula for evaluating the GMQD for any two-qubit 
state in 2010. Subsequently, Luo and Fu |l8[ introduced another form of GMQD and derived an explicit formula for 
evaluating the GMQD for an arbitrary state in a bipartite quantum system. However, using this way, usually, it is 
also difficult to calculate GMQD, and its explicit analytical expression for any bipartite system was not given. 

Besides the characterization and quantification of quantum correlations, another important problem is the behavior 
of these correlations under the action of decoherence. Recent results of the dynamics of the quantum correlation for 
a certain class of states in a two-qubit system under the influence of common noise channels show that the quantum 
correlations vanishes at asymptotic time. It maybe include a peculiar sudden change in behavior and is more resistent 
to the action of the environment than entanglement [l9|-[24j ■ Quantum correlations for multi- valued quantum system 
is a new and immature research area. Ali [1 61 ] studied the quantum discord for a two-parameter class of states in a 
2 <g> d system in 2010. For a qubit-qutrit (2 <g) 3) system, Ann [25| and Kan [26| studied its entanglement dynamics 
under the influence of dephasing and depolarizing channels, respectively. Ali et al. [27], HH and Ramzan et al. [29| 
studied its entanglement dynamics under phase damping and amplitude damping channels. Karpat et al. investigated 
its correlation dynamics under dephasing in 2011 }30( . All the states concerned in previous works are some classes of 
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the states in a qubit-qutrit system. 

In this paper, we devote to investigate the GMQD for a two-parameter class of states in a qubit-qutrit system and 
the dynamics of GMQD under the influence of various dissipative channels [i.e., both two independent (multi- local) 
and only one (local) dephasing, phase-flip, bit- (trit-) flip, bit- (trit-) phase-flip, and depolarizing channels] in the 
first time. Analytical results are presented. Our results show that all these dynamic evolutions do not lead to a 
sudden vanishing of GMQD. Quantum correlations vanish at an asymptotic time for local or multi-local dephasing, 
phase-flip, and depolarizing noise channels. However, it does not disappear even though t — > oo for local trit-flip and 
local trit-phase-flip channels. 



2. Initial states, noise model and geometric measure of quantum discord 

The class of states with real parameters in a hybrid qubit-qutrit (2 ® 3) quantum system [3l| are given by 

p bc (0) = a (|O2)(O2| + |12)(12|) + 6(|0+)(0+| + |0-)(0-| + |V^)^ + |)+c|^">(V^| ! (1) 

where 

I0±) = -L(|oo>±|ii», 

|^) = -^(|01) ±|10)), (2) 

and a, 6, and c are three real parameters, and they satisfy the relation 2a + 36 + c = 1. |0) and |1) are the two 
eigenstates of a two-level quantum system (qubit) or the eigenstates of a three-level quantum system (qutrit) with 
the other eigenstate |2). The two-parameter class of states pb c (0) can be obtained from an arbitrary state of a 2 <g) 3 
quantum system by means of local quantum operations and classical communication (LOCC) [31j ]. 

Let us briefly recall the quantum discord of a bipartite state p in a Hilbert space H A <g> H B . The quantum discord 
V(p AB ), a measures of all quantum correlations, is defined as the difference between the total correlation and the 
classical correlation [HI, that is, 

V(p AB ) = I(p AB )-C(p AB ), (3) 

where 

I(p AB ) = S(p A ) + S(p B ) - S(p AB ). (4) 

Here, I(p AB ) represents the quantum mutual information (the total amount of correlations) of the two-subsystem 
state p AB . p A ( B ^ = tr B (A)(p AB ) is the reduced density matrix for the subsystem A(B). S(p) = — tr(plog 2 p) is the 
von Neumann entropy of the system in the state p. The other quantity, C(p AB ), is interpreted [THHH as a measure 
of the classical correlation of the two subsystems AB in the state p AB and it is defined as the maximal information 
that one can obtain, for example, about B by performing the complete measurement Ilk on H A , 



C B {p AB ) = sup js(p B ) - £p fc S(p B l fe ) j , (5) 

where p s ' fe = -p^O^k <& lB)p AB (J^k <£> Ib) is the postmeasurement state of B after obtaining the outcome k on A with 
the probability Pk = ^((n^ ® Is)p AB (Hk <8> Ib))- is a set of one-dimensional projectors on H A . 

We note that for a general mixed state, the (one-side) classical correlation of Eq.© may be asymmetric with respect 
to the choice of subsystem to be measured (Ca(p) Cb(p)), that is, the quantum discord in Eq.Q is not symmetric 
(Da{p) 7^ Db{p))- However, it is known that Da(p),Db(p) > 0, and Da(p) = Db{p) — if and only if p is a 
classical-quantum state, and p has the following expression: 

P AB = Y. P ^)a(*\®\j)b(jI (6) 

i,j 

where P,j is a probability distribution, and \i) a and \j)b are the orthonormal bases of system A and B, respectively. 
GMQD: the geometric measure of quantum discord of the state p is given by [17| 

D G (p) = mm\\p- X \\ 2 , (7) 

x 
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where the minimum is taken over the set of zero-discord states [i.e., D{\) = 0] and \\p — xl| 2 = t r (p—x) 2 1S the square 
norm in the Hilbert-Schmidt space. 
For any two-qubit state 



P AB 



^[I®I + J2(x i a l ®I + y i I®<j l ) + *%2(rijcri®o-jj], (8) 

i=l i,j=l 

its GMQD is given by 

D G (p) = ^(\\X\\ 2 + \\R\\ 2 ~k max ), (9) 

where <7j are the Pauli spin matrices, X — (x\, X2, x^) T , R is the matrix with elements r^, and k max is the maximal 
eigenvalue of the matrix K = XX T + RR T . 

For a general bipartite system H A ® H B , if we choose orthonormal basis sets in local Hilbert-Schmidt spaces of 
Hcrmitian operators {X{\ and {Yj} (i = 1, • • ■ ,d\ and j = 1, • • ■ , d B ), any bipartite state can be written as 



= Y J CijX i ^Y j , (10) 



where c^- = tr(pXi ®Yj). Its GMQD can be expressed as [18 1 



£> G (p) = trlCC 1 ) - maxirMCC* J A J ), (11) 

A 

where C = (cjj) and the maximum is taken over all d,A X d^-dimensional matrices A = (cifci)- Qfci ~ ^(|fc)(fc|Xi). 
{|fc)} is any orthonormal base for H A . Here fc = 1, 2, • • • , d,A and i = 1, 2, • • • , d\. 

In our physical model of noise for a qubit-qutrit system (composed of a two-level subsystem A and a three-level 
subsystem B), the two subsystems interact with their environments independently. The evolved states of the initial 
density matrix of such a system when it is influenced by multi-local environments can be given compactly by 

i=i j=i 

Here, the operators E A and are the Kraus operators which are used to describe the noise channels acting on the 

qubit A and the qutrit B, respectively. They satisfy the completeness relations J2i ^i^Ef — I and F^F^ = I 
for all t. 



3. Geometric measure of quantum discord for qubit-qutrit systems 

It is important to consider the behavior of quantum correlations under the action of decoherence. In this section, 
we investigate what happens to the quantum correlations (i.e., GMQD) in a qubit-qutrit system under common noise 
channels for qubit (qutrit): dephasing, phase-flip, bit-(trit-) flip, bit-(trit-) phase-flip, and depolarizing channels. The 
Kraus operators describing these channels for a single qubit A and a single qutrit B are presented in Appendix. The 
time-dependent parameters are defined as 7^4 = 1 — e~ tTA and 7s = 1 — e~ trB , 7^4, jb & [0, 1] and Ta (Tb) denotes 
the decay rate of the subsystem A (B). The two specific environment noise situations will be considered: (i) local 
and (ii) multi-local. In the case (i), only one part of a qubit-qutrit system (S) interacts with its environment. In the 
case (ii), both the two parts of S interact with their local environments, independently. 

We choose 4 Hermitian matrices for each H A , 

Xl = 7f(ii)' X2 = 7i(?J)' X3 = 7l0?)' X4 = 7f(o-°0' (13) 
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and 9 matrices Hermitian matrices for each H B , 
Yi = 
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^(fl oj' 5 ^ \ 1 J ' 6 ^ 

! [1 \ /o \ j 

F 7 = — [ 1 I , Y 8 = — [ 1 I , Y 9 = —\ -/ I . (1.4) 

Any orthogonal base for H can written as 

|0„) = cos6>|0) +e^sin6»|l), 

\9_l) =sin^|O)-e^cos0|l). (15) 

The matrix A = (aki) in Eq. (|ll[) with the matrix elements: aki — tr(\k){k\Xi). Here, {|fc)} is any orthonormal base 
for H A shown in Eg. ([151) , an d Xi is a set of Hermitian matrices for H A shown in Eq. (fl~3l) . The matrix can be presented 
as 



. _ 1 / 1 sin 29 cos <p sin 29 sin <p cos 29 
y/2 \ 1 — sin 29 cos ip — sin 29 sin y> — cos 29 



(16) 



C-ZJ Channel without noise. The coefficient elements = tr(pX{ <8>Yj) of the matrix C in Eq. flTTj) are given by 

cn = 7!' 

2 - 96 - 3c 

Cl7 = — iTT - ' 

C22 = c 33 = c 44 = -(6- c), (17) 

and all the remaining matrix elements are zero. 

By replacing the factors A and C = (cy) in Eq. lTTl) with Eqs. dTB)) and ((TT]) . respectively (that is, Dq is the function 
of 9 and ip), and calculating the minimum of Dq over 9 and we obtain the GMQD for the qubit-qutrit systems 
under a channel without noise, that is 

D G {p AB ) = \{b-cf. (18) 
(2) Multi-local dephasing channel. The coefficient matrix elements cy are given by 



Cll 



1 

1 



ci7 = -^=(2 -96 -3c), 
C22 = c 33 = -(6 - c) v /(l - 7^)(1 - 7s), 

c 44 = ^(&-c), (19) 

and all the remaining matrix elements are zero. 

By replacing the factors A and C = (cy) in Eq. (fTTj) with Eqs. (fT6| and (fl9| . respectively, and calculating the 
minimum of £>g over 9 and <y9, we obtain the GMQD for the qubit-qutrit systems under a multi-local dephasing 
channel, that is 

D G { P AB )i = \{b - c) 2 (l - 7 a)(1 - 7b). (20) 
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Its dynamics is shown in Fig[T] In the paper, we consider the parameters c ^ b, that is, the initial state is a quantum 
correlation state. 

(3) Multi-local phase- flip channel. The coefficient matrix elements are given by 



en 



1 



ci7 = -^=(2 -96- 3c), 

C22 = c 33 = ^(6 - c)(l - 7 j4 )(1 - 7b), 

C44 = (21) 

and all the remaining matrix elements are zero. 

By replacing the factors A and C — (cy) in Eq. ([TT|) with Eqs. ([T6l) and (f2Tj) . respectively, and calculating the 
minimum of Dq over 9 and ip, we obtain the GMQD for the system, that is 

D G (p AB ) 2 = \{b - c) 2 (l - 1A f(l lB f. (22) 

Its dynamics is shown in Fig|2] 

(4) Multi-local bit- (trit-) flip channel. The coefficient matrix elements Cy are given by 



1 

71' 



ci7 = ^(2-96-3c)( 7 b-1) ! 

c 22 = -^(6-c)(2 7B -3), 

C25 = c 28 = ^(6- c)tb, 

c 33 = ^(6-c)(2 7fl -3)(7 A -l), 

C36 = -C 39 = ^(b-c)(j A -^hB, 


c 44 = i(6- C )( 7B -l)(7A-l), (23) 




FIG. 2: Dynamics of GMQD for the system undergoing the multi-local phase-flip and depolarizing noises. 



and all the remaining matrix elements are zero. 

By replacing the factors A and C = (cy) in Eq. (fTTj) with Eqs. ([T6|) and (f23|) . respectively, and calculating 
minimum of Dq over 9 and ip, we obtain the GMQD for the system, that is 

D G ( P AB h = ^(b - c) 2 (l - 7A ) 2 (6 + 5( 7B - 2) 7B ). 

Its dynamics is shown in Fig 13] 




FIG. 3: Dynamics of GMQD for the system undergoing the multi-local bit- (trit-) flip noise. 
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(5) Multi-local bit- (trit-) phase- flip channel. The coefficient matrix elements c^- are given by 



en 

C17 

C22 
C25 
C33 
C36 
C44 



1 

1 

V3 



(2-96-3c)( 7B - 1), 



-(6- c )(2 7s -3)( 7A -l), 

C28 = -^{b-c)(j A - l)7s, 

-i(&-c)(2 7B -3), 

-C39 = To( &_ C )7B, 
i(6-c)( 7B -l)( 7A -l), 



(25) 



and all the remaining matrix elements are zero. 

By replacing the factors A and C — (cy) in Eq. ([TTj) with Eqs. (IT6l) and (|25|) . respectively, and calculating the 
minimum of Dq over 9 and ip, we obtain the GMQD for the system, that is 



D G (p AB h = ^(b-rril - - ^-a-2 + lD m L! - 21)11. 



(26) 



Its dynamics is shown in Fig J?) 




FIG. 4: Dynamics of GMQD for the system undergoing the multi-local bit- (trit-) phase-flip noise. 



(6) Multi-local depolarizing channel. The coefficient matrix elements Cij are given by 

ci7 = ^=(2-%-3c)( 7s -l), 

C22 = c 33 = c 44 = ~(6 - c)(l - 7 a)(1 - 7 s). (27) 

all the remaining matrix elements are zero. 

By replacing the factors A and C = (c^) in Eq. ([PTj) with Eas. (IT6l) and (l27l) . respectively, and calculating the 
minimum of Dq over 9 and y>, we obtain the GMQD for the system, that is 

D G ( P AB h = \{b c) 2 (l - X4) 2 (l - IB?- (28) 



8 

Its dynamics is shown in Figj2] 

(7) Local qubit noise only. Consider 7^ = 0, the GMQD can be calculated as 

D£\p A % = \{b-cf{l- 1A)l 

D%\p AB ) 6 = \{b-cf{l- lA f. (29) 

Here D^q (p AB )§ corresponds to dephasing channel, and Dq (p AB )e corresponds to phase-flip, bit-flip, bit-phase-flip 
or depolarizing channels. The dynamics of GMQD for these cases are shown in FigJS] 




0.0 0.2 0.4 0.6 0.8 1.0 

7A 

FIG. 5: Dynamics of GMQD for the system undergoing the various local noises which act on the qubit alone. The solid and 
short-dashed lines correspond to dephasing and phase-flip (or bit-flip, bit-phase-flip, depolarizing) noises, respectively. 



(8) Local qutrit noise only. Consider 7^4 = 0, the GMQD can be calculated as 





2 V 


c) 2 (1-7b), 




D { S\P AB )7 = 


2 V 


C ) 2 (1-7b) 2 , 




D { *\p AB h = 


12 V 


-c) 2 (6 + 5( 7B -2)7 B ), 




Df{ P AB h = 


2V 6 


-c) 2 (12 + 75(975 -20)). 


(30) 



Here Dq\p ab )j, D^q [p AB )i, Dq\p ab )t, and Dq\p ab )7 corresponds to dephasing, phase-flip (or depolarizing), 
trit-flip, and trit-phase-flip channels, respectively. The dynamics of GMQD for these cases are shown in Fig|6l 

Eqs. ([29|) and (f30|) show that the environment, which acts on the qubit alone and causes phase- flip, bit-flip, bit-phase- 
flip, and depolarizing, affects the GMQD of a hybrid qubit-qutrit system in the same way. However, if the environment 
acts on the qutrit alone, only the phase-flip and depolarizing channels affect the GMQD of the qubit-qutrit system in 
the same way. 

4. Discussion and Conclusion 

Quantum discord T> for a two-parameter class of states in a hybrid qubit-qutrit system was discussed by Ali (l6j 
in 2010. We have studied GMQD and its dynamics under the influence of various dissipative channels, including 
local and multi-local dephasing, phase-flip, bit- (trit-) flip, bit- (trit-) phase-flip, and depolarizing channels, in the 
first time. Moreover, the explicit analytical expressions were gave out. Our results show that environment, which 
causes dephasing, phase-flip, bit- (trit-) flip, bit- (trit-) phase-flip, and depolarizing of a qubit-qutrit system, affects 
the quantum correlations of a hybrid qubit-qutrit system in very different ways. All these dynamic evolutions do 



9 




FIG. 6: Dynamics of GMQD for the system undergoing the various local noises which act on the qutrit alone. The solid, dashed, 
short-dashed and dashed-dotted lines correspond to dephasing, phase-flip ( or depolarizing), trit-flip, and trit- phase-flip noises, 
respectively. 



not lead to a sudden vanishing of GMQD. Quantum correlations vanish at asymptotic time for local or multi- local 
dephasing, phase-flip, and depolarizing noise channels, while it cannot be destroyed completely even though t — > oo 
for local trit-flip and local trit-phase-flip channels. The states shown in Eq.([T]) with a = (that is, c = 1 — 36) are 
equivalent to Werner states 35] in 2 <g> 2 system with different noise channels. Compared with the states in two-qubit 
systems, their quantum correlation all vanishes at a asymptotic time and can not occurs quantum correlation sudden 
death and sudden birth under an arbitrary Markovian dynamics [l9l [2ll [22| . This phenomenon maybe provide some 
important information for the application of GMQD in hybrid qubit-qutrit systems in quantum information. 
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Appendix A: Various dissipative channels 



Dephasing channels: the set of Kraus operators for a single qubit A that reproduces the effect of a dephasing 
channel are given by [33[ 



*A I i o \ T „ A I 



^ = ( o VT^Xr h > E * 1° Vtz'® 73 ' (A1) 



and those for a single qutrit B can be written as [25j 

1 \ / 

F? =h®\ yr^i , F 2 B = I 2 ® ) , F 3 B =/ 2 ®| ] . (A2) 

V1-7B / I 



The time-dependent parameters are denned as 7^4 = 1 — e~ tTA and 7s = 1 — e~ tTn . Here 7^4, 75 £ [0, 1]. Fa (Tb) 
denotes the decay rate of the subsystem A (B). 

Phase-flip channels: the Kraus operators describing the phase-flip channel for a single qubit A are given by [33| 




sf^/i-f J? ® J * # = Vf U-°i)® J »' (A3) 
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and those for a single qutrit B can be written as 



1 

e ~ l27r / 3 
e' i27r / 3 



1 

e i27r / 3 
e -^/3 



> (A4) 



where 7^ = 1 — e tVA , 75 = 1 — e ITb , and 7a, 7s G [0, 1]. (T^) represents the decay rate of the subsystem A 

(B)- 

Bit-(Trit-) flip channels: the Kraus operators describing the bit- flip channel for a single qubit A are given by [33J 



£f = Wl- 



7a ( 1 
1 



E A = 



7 A I 1 
1 



and those for a single qutrit B can be written as 

Ff = 
where 7^4 = 1 






(A5) 



(A6) 



-tTi 



, 7b = l-e~ tlB , and 7,4,75 € [0,1]. 



Bit-(Trit-) phase-Hip channels: the Kraus operators describing the bit-phase flip channel for a single qubit A are 
given by 33] 



pA _ L 1A ( 1 

Si - V Tlo 1 

and those for a single qutrit -B can be written as 



F" = h 



2 7i 



F 2 B = h 



7A ( Q —i 



i 



(A7) 



e -^/3 

1 
e i27r/3 

e l27r / 3 

e - i2 V3 

1 



F 



e i27r / 3 
10 
e ~ l2jr / 3 

e - i27r / 3 

e j27r / 3 

1 



(A8) 



where 7^4 = 1 — e A , 75 = 1 



,-tr* 



and 7a,7b G [0, !]• 



Depolarizing channels: the set of Kraus operators that reproduces the effect of the depolarizing channel for a single 
qubit A are given by [H[ 



(A9) 



where <r, (i = 1,2,3) are the three Pauli matrices. The Kraus operators describing a single-qutrit depolarizing noise 
are given by [? ] 



F\ 



F? 



Fy = 



h 

h 



V 1 -— /3 > 


Fi 


= h <X 


3 ' 


Fi = 


h ® 


3 ' 
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Fi 


= J 2 ® 




Fi = 




^V 2 z, 


3 


Fi 


= I 2 <S> 


3 


Fi = 


= 7 2 (8 


,V^2 

3 



where 



Y 





1 






• 
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■-( 












1 

e t2rr / 3 







-i2w/3 



(A10) 



(All) 
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and 7a = 1 - e~ trA , 7s = 1 - e - tr «, 7a , 7b e [0, 1]. 
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